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ABSTRACT 

We present a measurement of the surface brightness of the cosmic X-ray background (CXB) in the Chandra 
Deep Fields, after excluding all detected X-ray, optical and infrared sources. The work is motivated by a recent 
X-ray stacking analysis by Worsley and collaborators, which showed that galaxies detected by HSThvX not by 
Chandra may account for most of the unresolved CXB at £ > 1 keV. We find that after excluding HST and 
Spitzer IRAC sources, some CXB still remains, but it is marginally significant: (3.4 ± 1.4) x 10"'-^ ergs cm"^ 
s"' deg"2 in the 1-2 keV band and (4 ± 9) x 10"'^ ergs cm"^ s"' deg"^ in the 2-5 keV band, or 7% ± 3% and 
4% ± 9% of the total CXB, respectively. Of the 1-2 keV signal resolved by the HST sources, 34% ±2% comes 
from objects with optical colors typical of "normal" galaxies (which make up 25% of the HST sources), while 
the remaining flux comes from objects with colors of starburst and irregular galaxies. In the 0.65-1 keV band 
(just above the bright Galactic O VII line) the remaining diffuse intensity is (1.0 ±0.2) x 10"^^ ergs cm~^ s~^ 
deg~^. This flux includes emission from the Galaxy as well as from the hypothetical warm-hot intergalactic 
medium (WHIM), and provides a conservative upper limit on the WHIM signal that comes interestingly close 
to theoretical predictions. 

Subject headings: methods: data analysis — X-rays: diffuse background — X-rays: galaxies 



1. INTRODUCTION 

The cosmic X-ray background (CXB) is known to be pri- 
marily the integrated emission from X-ra y point sources , 
mostly active ga lactic nuclei (e.g., Hasinger et alj 120051 : 
iBauer et ani2004l) . Recent studies have shown that the Chan- 
dra Deep Fields (CDFs) North and South, the most sensi- 
tive observations of the X-ray sky to date, resolve - ^80% 
of die extr agalactic CX B at 1-2 k eV (Moretti et all 120031: 
I Worsley et al. 2005; Hi ckox & Markevitch 2006i hereafter 
HM06), and -50% at 7 keV (IWorsley et al.ll2005h . While 
the majority of the CXB at £ < 5 keV has been accounted 
for, the nature of the still-unreso lved com ponent (the absolute 
intensity of which we derived in lHM06h remains unclear. 

Part of the unresolved CXB is due to faint p oint sour ces that 
have yet to be detected in X-rays. However. IHM06| showed 
that an extrapolation of the observed logA^/ log5 curve to zero 
fluxes falls far short of the unresolved CXB. To account for it, 
the slope of the relation must steepen below the CDF-N X-ray 
source detection limit of 5'o.5-2keV = 2 . 4 x 10 "'^ ergs cm"^ s"' . 
This is not implausible; iBauer et aP (l2004l) showed that the 
CDFs contain a population of normal and starburst galaxies 
with 5'().5-2keV < 10"'^ ergs cm"^ s"' whose logA^/logS' is in- 
deed steeper, and which should dominate the source counts at 
So.5-2 keV < 10"'^ ergs cm~^ s~' . These faint sources, while not 
detected in X-rays, can be associated with optical or infrared 
sources found in the CDFs. A recent X-ray stacking analy- 
sis of Hubble Space Telescope (HST) sources in the CDFs has 
found that their X-ray emission can account for most of the 
unresolved E = 1 - 6 keV CXB (Worsley et al. 2006). 

Still, some of the unresolved flux should be truly dif- 
fuse. At £ < 1 ke V, there is flux from the Local Bub- 
ble dSnowdenl |2004|) as we l l as charge exchang e local to 
the Sun (e.g . faavensi 120001: IWargelin et al.ll2004l ). In addi- 
tion, some signal may come from the warm-hot intergalac- 
tic medium (WHIM) with temperatures 10^-10^ K, that is 
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thought to cor nprise the m ajority of the baryons in th e local 
Universe (e.g., ICen feOstr iker 1999i: lDave etalj2001h . Mod- 
els of the WHIM predict a signal dominated by line emission 
from oxygen and iron with p rominent lines at . 5-0.8 keV, 
but some fl ux up to 1.5 keV (iRiimps et al.ll2001l: iFang et alj 
12005; Urs ino & Galeazzil l20(m To date the WHIM has 
proved challenging to observe, although there have been re- 
ports of detections of WHIM X-ray absorption lines (e.g, 
iNicastro et al.ll200"5l however see Kaastra et al. 2006), as well 
as excess soft X-ray emission around galaxies that has been 
attributed to the WHIM (e.g., Soltan 2006). Diffuse emission 
could also come from mo re exotic sources such as decay of 
dark matter particles (e.g., lAbazaiian et"al]|2007h . 

In this Letter, we constrain the diffuse CXB by directly 
measuring its absolute intensity in deep Chandra pointings 
after the removal of all sources detected in deep X-ray, opti- 
cal, an d IR observations. Our analysis is complementary to 
that of IWorsley et all (12006) . and is a foflow-up to HM06u in 
which we excluded only the X-ray sources. 

2. data 

Chandra X-ray observations of the CDFs North and South 
(CDF-N and CDF-S) have respective total exposure times of 
roughly 2 Ms and 1 Ms. The CDF-N data consist of two sub- 
sets, observed in in Very Faint (VF) or Faint (F) telemetry 
mode. Each of the three subsets (CDF-N VF, CDF-N F and 
CDF-S) contains 6 to 9 individual observations with exposure 
times 3 to 17 ks each. The datasets are identical to those 
used in lHM06l except for the ACIS calibration updates. For 
subtraction of the instrumental background, we add the latest 
ACIS stowed observations, for a total background exposure 
of 325 ks, compared to 256 ks used in HM06. These updates 
caused no significant change to the CXB fluxes. 

In HM06, we measured the CXB intensity after excluding 
X-ray detected sources from the catalog of Alexander et alj 
(2003). To further exclude the contributions of X-ray faint 
objects, we use sources detected in the optical and IR ob- 
servations of the Great Observa tories Origins Deep Survey 
(GOODS, iDickinson etani2003h . Optical and near-IR data 
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Fig. 1 . — Smoothed 0.5-8 keV image of CDF-N, showing the source ex- 
clusion regions. The large circle is the 3.2' extraction region, medium-size 
circles and ellipses show X-ray point and extended sources, and small circles 
show HST sources. 
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Fig. 2. — CXB spectrum for the CDF-N VF subset, with all sources (Chan- 
dra, HST and IRAC) excluded. The spectrum is fit with an APEC model with 
kT = 0.20 keV plus a power law with F = 1.5. Black points show the data in- 
cluded in the fit (se dHMOffl , and the red dashed line shows the best -fit model 
for the CDF-N VF spectrum excluding only X-ray sources )HM06t) . 



were taken using the Advanced Camera fo r Surveys (ACS) 
on HST in the B435, Veoe, '775, and zs50 bands CiGiavalisco et alj 
|2004). We use the pubHc catalog^ based on detections in the 
Z850 band (A ^ 8300-9500 A), with a Hmiting magnitude of 
Z850 ^ 27. The //5r catalogs for CDF-N and CDF-S each con- 
tain ~4500 objects within the area for which we extract CXB 
spectra (inside a 3.2' circle, and away from detected X-ray 
sources, see §|3]l. 

We also exclude sources detected with the Infrared Array 
Camera (IRAC) on Spitzer, in four bands centered on 3.6, 4.5, 
5.8, and 8 /im. There is no published catalog currently avail- 
able for these sources, so we have performed our own source 
detection using the calibrated IRAC imag es-^. Sources were 
detected using the SExtractor code dBertin & ArnoutsI 
1 19961) in all four bands. In the most sensitive band (3.6 //m), 
we detect ^1000 sources in the X-ray extraction region to 
a 5(7 Hmiting flux of -0.2 /iJy. Most (>95%) of the IRAC 
sources have counterparts within 1" in the HST zm) catalog. 

We also considered Spitzer Multiband Imaging Photometer 
(MIPS) 24 /im sourc es^ and for CD F-N, Very Large Array 
(VLA) 21 cm sources (lRichardsl2000l) . with limiting fluxes of 
80 and 40 /iJy, respectively. All the sources within our X-ray 
extraction regions have HST or IRAC counterparts, except for 
2^ of the 80-90 MIPS sources and 1 of the 6 VLA sources, 
so we do not include these data in the analysis. 

3. ANALYSIS 

We calcul ate the C XB surface brightness as discussed in de- 
tail in § 5 of HM06, with only minor differences. For the HST 
and IRAC sources not detected in X-rays, we use smaller ex- 
clusion regions with radii equal to rgo, the X-ray point-spread 
function (PSF) 90% enclosed flux radius at 1 .5 ke V (see Eqn. 
1 of HM06), compared to 4.5-9r9o for X-ray detected sources. 
Due to the large number of optical and IR sources, at large off- 
axis angles there is little area remaining after source exclu- 
sion. Therefore we only extract spectra from a circle around 
the optical axis with radius 3.2' (where rgo = 2.2"), compared 
to 5' in,HM06 . The solid angle we use here is 0.45 of that in 
lHM06l when we exclude only X-ray sources, and 0.32 when 

^ http://archive.stsci.edu/pub/hlsp/goods/catalog_rl/h_rl.lz_readme.html 
' http://data.spitzer.caltech.edu/popular/goods/Documents 
/goods_dataproducts.html 



TABLE 1 

CXB SURFACE BRIGHTNESS AFTER EXCLUDING ALL SOURCES 



Energy 




CXB intensity" 




CXB frac.'' 


(keV) 


CDF-S 


CDF-N (VF-l-F) 


Average 


(%) 


0.5-2=^ 


45 ±8 


41±11'' 


43 ±7'' 




0.65-1'-' 


11±2 




10 ±2'* 




1-2 


5.0±1.8 


2.4±1.6 


3.4±1.4 


7.3 ±3.0 


1-1.3 


1.7±0.7 


1.0±0.7 


1.3±0.6 


8.3 ±3.8 


2-5 


7±12 


2 ±10 


4±9 


4±9 


2-8 


14 ±39 


4±30 


7±26 


4±16 


" Units are 




;s cm - s deg , 


errors are 68 


%.•> Percent 



of total CXB, assuming a power-law CXB spectrum with F = 1.4 
and nor malization 10.9 ±0.5 photons cm"^ s"' sr"' keV"' at 1 keV 
<HM0d) .'^ Intensities for E < I keV are dominated by Galactic and 
local line emission that varies with time and position. These are not 
extragalactic values. Does not include CDF-N VF, to exclude the 
local temporary excess. 

we exclude all HST and IRAC sources. 

For each of the datasets, we fit the spectrum of the remain- 
ing area with a model consisting of an optically-thin thermal 
component (APEC) with solar abundances and kT ~ 0. 15-0.2 
keV, to account for emission from the Galaxy as well as any 
diffuse thermal CXB. We also include a power law compo- 
nent that dominates at £ > 1 ke V, an d fix the photon index at 
r = 1.5 , the best-fit value from HM06I (we do not vary F as in 
IHM061 due to poorer statistics). These models do not neces- 
sarily provide a physical description of the emission, but they 
are sufficient for our purpose of calculating fluxes. A CDF-N 
VF spectrum with the model fit is shown in Fig.|2l 

For each spectrum, we calculate fluxes in the 0.5-2 keV, 
1-2 keV, 2-8 keV, and 2-5 keV bands, as well as the 0.65- 
1 keV and 1-1.3 keV bands, where the hypothetical WHIM 
signal or galaxy group emission are be st di stinguished from 
the diffuse Galactic foreground (see § I5.2l i. Errors on the 
observed flux include statistical errors in the sky and back- 
ground count rates, and a 2% sys tematic uncertainty in the 
background normalization (§ 4 of lHM06l) . We carefully ac- 
count for background uncertainties. Because of differences 
in pointing of up to 2' and the small X-ray extraction area, 
the composite stowed background spectrum includes 3 times 
(for CDF-N) and 1.5 times (for CDF-S) as many statistically 
independent photons as do the spectra for each separate ob- 
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Fig. 3. — (a) Solid angles of the X-ray extraction regions, and (b-d) unresolved CXB surface brightness, for subsets of excluded Chandra, HST, and IRAC 
sources. In all panels, points on the far left show the value when only X-ray sources are excluded (as in HM06, but using the smaller 3.2 ' region). Points with 
Z850 labels on the x-axis correspond to the additional exclusion of HST sources brighter than z.s5o (note that z.g5o = 27 is only an approximate detection limit for all 
HST sources). Points on the right show intensities excluding X-ray sources plus only //ST "normal" galaxies (§ 15. IK only IRAC 3.6 ^m sources, and finally all 
HST and IRAC sources. For the 0.65-1 keV band we show the CDF-N VF and CDF-N F subsets separately, as time variability in the local emission (apparently 
charge exchange in the Solar wind, e.g., Wargelin et al. 2004) caused these intensities to differ. For clarity, errors are shown only on the far left and right points. 
For (c) and (d), on the right axis we show the corresponding fraction of the total CXB, for which we assume a power law spectrum with F = 1 .4 and normalization 
10.9 photons cmr^ s"' sr"' at 1 keV jHM06l) . 



servation. We accordingly increase the effective background 
exposure times by these factors. We conservatively neglected 
this effect in HM06, as it was much smaller due to the larger 
extraction area and fewer excluded sources. 

A further small correction is necessary to subtract the point 
source flux that is scattered outside our exclusion regions. For 
each HST or IRAC source, 10% of the X-ray flux should lie 
outside rgQ. However, because many regions overlap, the con- 
tribution of this scattered flux is smaller Using a model of 
the Chandra PSF, we estimate that when all HST and IRAC 
sources are excluded with radius ryo, 5% (0.5-2 keV) and 8% 
(2-8 keV) of their flux should remain in our unresolved spec- 
trum. Likewise, 0.3% and 0.6% of the flux from detected X- 
ray sources (for which we use larger exclusion regions) will be 
scattered into the spectrum. We account for both these small 
contributions in the final unresolved intensities. As a check, 
we tried exclusion regions of radius 1 .Srgo and obtained iden- 
tical results. 

4. RESULTS 

We first calculate the CXB intensity after excluding HST 
subsamples with various limiting magnitudes from zaso < 21 - 
26, and finally excluding all HST sources (z < 27). We also 
try excluding only "normal" galaxies, selected based on their 
HST colors (see § 15. Il l, only IRAC 3.6 /im sources, and finally 
all HST and IRAC sources. In every case we also exclude all 
detected X-ray sources. The results are shown in Fig.|3] and 
the residual CXB intensities excluding all sources are given 



in Table 1 . We also calculate the mean CXB intensity for all 
three subset s, includ ing a detailed error analysis as described 
in § 8.2.3 of lHM06l 

Fig. [3] shows that excluding faint HST sources has little 
effect on the CXB at £ < 1 keV, where it is dominated by 
the truly diffuse Galactic and local emission. At £ > 1 keV, 
as fainter objects are excluded, the surface brightness of the 
CXB decreases, indicating that the HST sources do indeed 
contribute to the X-ray unresolved CXB. In the 1-2 keV band, 
the remaining intensity after excluding all HST and IRAC 
sources is above zero at only 1.5ct for CDF-N and 2.8(t for 
CDF-S. The CXB intensity decreases even as we excluded 
the faintest HST sources, suggesting that if we were able to 
go slightly deeper, we would resolve the 1-2 keV CXB to 
within our absolute uncertainties. The X-ray undetected HST 
sources contribute 1-2 keV intensities of (5.7 ±0.3) x 10"'^ 
and (4.3 ±0.3) x 10"'^ ergs cm'^ s"' deg-^ in CDF-N and 
CDF-S, respectively, or 59% ± 9% and 45% ± 7% of the un- 
resolved CXB from lHM06i The fluxes resolved by HST 
sources (given here and in § 15.11 ) represent differences be- 
tween the points in Fig. [3] and do not depend strongly on 
the instrumental background. Thus their errors are much 
smaller than the uncertainties on the absolute surface bright- 
ness shown i n Fig. [3] Stacking the X-rays from the same 
HST sources. IWorsley et alJ (l2006h found larger intensities of 
(8.4 ±0.5) X 10-1^ and (9.3 ±0.7) x 10"'^ ergs cm'^ s"' deg-^ 
for CDF-N and CDF-S, respectively. The cause of this differ- 
ence is not immediately clear, but may be related to different 
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Fig. 4. — Color selection of "normal" galaxies from HST data. All points 
show B435 — ^606 vs. !775 - ztsso colors for sources in the UDF with zs50 < 27. 
Sources are classified as elliptical and spiral galaxies (red circles) or starburst 
and Im galaxies (black dots) from photo-z analysis (Coe et al. 2006). Lines 
show our criteria to select "normal" galaxies from the GOODS data. 

accounting for overlap between the HST source regions. 

In the 2-5 keV band, the HST and IRAC sources account 
for the entire unresolv ed CXB within the large absolute un- 
certainties. Similarly, IWorslev et al.l (|2006) found that the 
stacked flux from HST sources is consistent with completely 
resolving the 2-6 keV CXB. 

5. DISCUSSION 
5.1. Which sources resolve the CXB? 

iBauer et al.1 (|2004 and lHM06l show that AGN likely con- 
tribute a relatively small fraction of the CXB at X-ray fluxes 
below the CDF limit, where the X-ray source population is 
dominated by starburst and normal galaxies. To determine 
what class of galaxies gives the largest contribution to the un- 
resolved CXB, we use HST colors to approximately divide 
the Z850 sources into two subsets by galaxy type. While we do 
not have optical classifications for the HST GOODS galax- 
ies, such a classification exists for the HST Ultra Deep Field 
(UDF), which covers a small region roughly 3' from the CDF- 
S aimpoint with a flux limit of z%sq — 30. The UDF was 
observed i n the same four ACS bands as the wider GOODS 
fields, and lCoe et al.l (l2006h performed a detailed photometric 
redshift analysis of the UDF sources, which classified each 
galaxy with one of eight templates (elliptical. Sab, Sbc, Im, 
and four types of starburst). 

In F ig, m we show B435 -V'606 vs. ins-zm) colors for 
ICoe et al. (2006) sources with mo < 27, corresponding to the 
flux limits of the larger GOODS fields. In this diagram, nor- 
mal (elliptical and spiral) galaxies occupy a different region 
from starburst and Im galaxies, with approximate boundaries 
'775 -Z850 = 0.2 and B435-V6O6 = 0.9- Oto -Z85o), as shown 
in Fig. m Of the 255 "normal" galaxies in the UDF sample, 
201 (73%) are in this region, while only 52 (21%) of the 253 
galaxies in this region have starburst or Im classifications and 
would be mis-classified by this criterion. In the GOODS cat- 
alog, these criteria classify 25% of zgso sources as "normal" 
galaxies and 70% as "starburst/Im" objects (5% of the sources 
do not have detections in all four ACS bands). The "normal" 
galaxies contribute 34% ± 2% of the 1-2 keV intensity re- 
solved by all the HST sources, and so have higher average 
X-ray fluxes than other X-ray undetected HST objects. 



We also examine the contribution from the 3.6 fim sources, 
which are on average only 0.25 magnitudes brighter in zggo 
than the HST sources, but account for a disproportionately 
large fraction of the unresolved CXB. The ~1000 3.6 /im 
sources contribute 63% ±2% of the total 1-2 keV CXB re- 
solved by all HST and IRAC sources. In comparison, the 1000 
brightest HST sources (which have zsso < 24.5, and of which 
only 20% are IRAC sources) account for 53% ±2%. There- 
fore, the X-ray flux from undetected X-ray sources is appar- 
ently more highly correlated with 3.6 /im flux than with zggo 
flux. This may be an important clue to the nature of these 
faint X-ray sources, although a detailed discussion is beyond 
the scope of this Letter. In a forthcoming paper, we will de- 
termine how bright these undetected sources can be in the X- 
rays and whether Chandra could detect them in a reasonable 
longer exposure. 

5.2. Comparison to WHIM predictions 

While most of the remaining CXB signal should be due 
to Galactic or local diffuse emission and unresolved X-ray 
point sources, our results provide conservative upper limits on 
emission from the WHIM (Cen et al. 2000 and later works). 
The brightest WHIM emission feature should be a redshifted 
O VII 0.57 keV line, but it cannot be distinguished from 
the bright Galactic O VII line without much better ('^l eV) 
spectral resolution. However, the hotter WHIM phases with 
kT ^ 0.2-0.5 keV would also emit significant Fe XVII lines 
at £ « 0.72-0.83 keV (rest-frame). Thus, the ratio of WHIM 
emission to Galactic foreground (which is well-approximated 
hy a kT ^ 0.15-0.2 keV thermal spectrum. Fig. 1) may be 
higher at £ > 0.6 keV. 

We can compare our 0.6 5-1 keV s ignal to predictions from 
cosmological simulations. iFang et a l. (2005) find WHIM in- 
tensities of ~ 3 X 10"**, 3 X 10"'", and 3 x lO"'"* ergs cm~^ s~' 
deg~^ toward a galaxy g roup, filament, and void, respectively. 
iRoncarelli et al.l (120061) obtained a sky-averaged 0.65-1 keV 
signal of 5 X 10"'^ ergs cm"^ s"' deg"^ after excluding galaxy 
groups and clusters. (For these two works, we have approx- 
imately converted fluxes from 0.1-1 keV and 0.5-2 keV to 
0.65-1 keV). The CDFs were selected to lie away from bright 
groups, but if they do not lie in voids, then our observed 
signal, (1.0 ±0.2) x 10"'^ ergs cm"^ s"' deg~^, comes inter- 
estingly close to the above prediction s. Galactic superwind s 
may increase the WHIM metallicity (ICen & Ostriker 2006*). 
so the line emission may be even higher. At higher energies, 
our 1-2 keV signal is 15 times higher than the predicted sky- 
averaged WHIM intensity of 0.19 x lO"'-' ergs cm"^ s"' deg"^ 
(Phillips et al. 2001), so does not give an interesting limit. For 
a more quantitative comparison, it will be necessary to excise 
the regions around individual galaxies in the WHIM simu- 
lations similarly to our analysis. With such a simulation, one 
may be able to place interesting constraints on the distribution 
and metallicity of the WHIM in the CDFs. 
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